Unambiguous determination of the g-factor for holes in bismuth at high B/T 
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Magnetotransport has been investigated in high purity bismuth crystals in static magnetic fields 
as high as 20 T and temperatures as low as 25 mK. This high B/T ratio permits observation of 
pronounced Shubnikov-de Haas oscillations over a wide field range and up to fields where most of the 
carriers are in the lowest Landau level. For transport currents in the bisectrix or binary directions 
and field in the perpendicular trigonal direction, we have observed doublet splittings centered on 
each Shubnikov-de Haas oscillation. These splittings exhibit a quadratic dependence on field and 
disappear before the last oscillation. Our observations allow us to conclude unambiguously that 
already when the Landau level index is as high as 2 the carriers are fully polarized and that the 
g- factor for holes with the field in the trigonal direction is 35.3(4). 



In this paper we report on magnetotransport measure- 
ments of crystalline bismuth at high magnetic fields and 
low temperatures. These measurements are motivated in 
part by the unusual electronic properties of the semimetal 
bismuth that reflect its unique location in an intermedi- 
ate position between good metals and semiconductors. 
The rhombohedral unit cell of bismuth can be obtained 
by slight deformation of a double cubic cell and con- 
tains two atoms with five electrons each. These ten elec- 
trons fill the first five Brillouin zones and then spill into 
the sixth zone leaving behind hole pockets in the fifth 
zone. The Fermi surface is well described by three ellip- 
soids containing electrons lying almost in the bisectrix- 
binary plane and one ellipsoid containing holes and ori- 
ented along the trigonal axis. The small cross section 
of this highly anisotropic Fermi surface leads to a small 
carrier density that is the same for both electrons and 
holes (n holc = n c icctron ~ 3 x 10 17 cm -3 ). Accordingly, 
bismuth is a perfectly compensated semimetal. This un- 
usual electronic structure together with the ready avail- 
ability of high quality crystals with lupg mean free paths 
has generated a large field of researchtn 2 ] on magnetic field 
induced oscillations in the magnetization (de Haas-van 
Alphcn effect) and conductivity (Shubnikov-de Haas ef- 
fect). 

For the highest magnetic fields, the carriers can be 
placed into the lowest Landau level and the quantum 
limit realized. The direction of the field is important. 
For example, if the field is along the bisectrix direc- 
tion (aligned closely with the electron ellipsoids) then 
the quantum limit can be realized at a few Tesla for 
the light electrons while the heavier holes remain in the 
quasiclassical regime in fields up to 50 T. Guided by 
these considerations and cognizant of the very limited 
experimental evidence suggesting the existence of high 
magnetic-field-induced correlated states,[J~El we have ini- 
tiated magnetotransport studies of bismuth at high B/T 
with the aim of finding magnetic field induced instabili- 
ties that might drive bismuth into a charge-density-wave 
(CDW) insulating state. Although our initial results do 
not show any evidence of such instabilities, they have 
clarified the meaning of doublet structures that coincide 



with the Shubnikov-de Haas oscillations and which have 
a heretofore- unobserved spacing that is quadratic in field. 
The observations and analysis reported here reveal an un- 
ambiguous determination of the g-factor or equivalcntly 
the ratio of the Zeeman energy to the cyclotron energy 
for holes confined to Landau orbits by field directed along 
the trigonal axis. 

The bismuth crystals used in this study were obtained 
from a variety of sources and had purities of 99.9995% or 
better. The crystals were cleaved at liquid nitrogen tem- 
peratures to expose the trigonal plane and then aligned 
along the bisectrix or binary directions according to 
the striation marks created during the low-temperature 
cleave. These alignments were confirmed using Lauc 
diffraction. Typical crystals were cut with a string saw 
to 5 mm length and had cross sectional areas on the or- 
der of 4 mm 2 . After cutting, the crystals were etched in 
a nitric/acetic acid solution and then annealed to tem- 
peratures as high as 250 C for periods as long as a few 
days. The small area contacts made to the samples using 
gold wires and low melting temperature solder (Woods 
metal) had contact resistances on the order of 0.1 Vt. The 
current direction was chosen to be along the binary (bi- 
sectrix) direction and the samples were mounted on a ro- 
tating stage so that the applied magnetic field could be 
rotated in the plane determined by the trigonal and bi- 
nary (bisectrix) directions. All measurements were taken 
using low frequency (< 20 Hz) ac techniques. 

Shown in Fig. (|l|) is the dependence of the longitudi- 
nal resistance on the reciprocal magnetic field for a binary 
sample with the magnetic field aligned perpendicular to 
the trigonal plane. The temperature is at 25 mK and 
the magnetic field was swept up to 20 T. This high B/T 
allows more than 30 Shubnikov-de Haas oscillations to 
be clearly discerned. For clarity, double logarithmic axes 
were used to accommodate the greater than two-decade 
range in field and greater than three-decade range in re- 
sistance. The choice of logarithmic axes however obscures 
the presence of almost uniform spacing A(l/B) of the 
Shubnikov-de Haas oscillations. 

Assigning consecutive integers n to the resistance min- 
ima of Fig. (j|) and then plotting the values of 1/B at the 
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FIG. 1. Resistance vs reciprocal field on logarithmic axes. 
The upper right inset is an expanded view of the n — 5 
through n — 10 oscillations and the lower left inset is an 
expanded view of the n = 1 oscillation showing the sharpness 
of the resistance change as the Fermi energy passes through 
the fully polarized v = 1 Landau level. 



respective minima versus n addresses this aspect of the 
data analysis. The plot of Fig. (|J) shows the results of 
such a plot for a set of integers that are consistent with 
a smooth extrapolation of the data to n = 0. The set 
of integers assigned in Fig. (fj) and labeled for n = 1 in 
Fig. (Fy) uniquely identifies the order of each Shubnikov- 
de Haas oscillation. As will be shown below, when spin 
splitting is taken into account, the order of oscillation n 
need not be the same as the the integer v identifying the 
order of the Landau level. 

At high fields the local slope in Fig. (Q) gives a pe- 
riod for the Shubnikov-de Haas oscillations A(1/_B) = 
0.146 T _1 , in good agreeme«t-|With experimental values 
obtained by others for holestffiJ with the field applied in 
the trigonal direction. The positive curvature reflects the 
fact that the valence and conduction bands are affected 
by the field, giving rise to a change in the carrier density 
and thus in the Fermi energy. For the field along the 
trigonal axis both the carrier density and the Fermi en- 
ergy increase with the fieldB This field-induced increase 
in the Fermi energy and the concomitant decrease in the 
1/B spacings of the Shubnikov-de Haas oscillations qual- 
itatively explains the positive curvature. We note that 
this effect is actually quite small when compared to what 
is expected when the field is applied in the basal plane, 
where the Fermi energy has a stronger dependence on B. 

A salient feature of the data presented in Fig. (|l|) is the 
appearance of doublets coinciding with the valleys of the 
Shubnikov-de Haas oscillations. These doublets persist 
to low fields as shown in the inset and as many as 10 
can be clearly identified. In addition the spacing in the 
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FIG. 2. Plot of the reciprocal fields as a function of integer 
n for the first twenty-five resistance minima of Fig. 1. The 
integers n identify the order of the Shubnikov oscillation with 
n = 1 labeled in Fig. 1 

field between the subpeaks of a doublet increases with the 
field. This observation is made more quantitative in the 
plot of Fig. (^) where the doublet spacings are plotted 
against the square of the field. 

For the oscillatory contribution to the holes' conduc- 
tivity for the field along the trigonal axis, one haalil 
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is the temperature smearing factor, and Rd{p) is the dis- 
order smearing factor. The parameters a, are the compo- 
nents of the inverse mass tensor with respect to the crys- 
tallographic axis (1 corresponds to the binary axis, 2 to 
the bisectrix one and 3 to the trigonal axes: for the holes, 
Qii = 0^2)7 mo is the bare electron mass, tp — 0.012 eV 
the hole Fermi energy at zero field, g the hole gyromag- 
netic ratio and Ib — yK/eB the magnetic length. The 
disorder smearing factor has the form 



Rd(T) = exp 
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where — h/2nT is the Dingle "temperature" and t 
is the hole mean free time. Using Eqn. (1) we analyzed 
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FIG. 3. Plot of the doublet spacings against the square of 
the field. The inset is a schematic of the energy levels showing 
the spin direction and the identifications (n, v) of the order 
of oscillation n and Landau level v. 



the smearing of the Shubnikov-de Haas oscillations at the 
relatively high temperature of 5K where the doublets are 
washed out and the first harmonic (p = 1) dominates. 
At this temperature all of the broadening is found to be 
thermal, implying that Trj is well below IK. This esti- 
mate is consistent with our experimentally determined 
estimate of r (t~2x 1CP 10 s) and corresponding lower 
bound on Td {Td > 0.006K) calculated using the Drudc 
formula and the measured resistivity at low temperature 
and zero field. 

Eqn. (jib, gives us the observed period with the mag- 
netic fieldE The signature of the crossing of a Landau 
level by the Fermi energy in a perfect crystal at zero tem- 
perature is a singularity in the oscillations and a sharp 
peak in a real crystal at high B/T ratio. In the ex- 
periment, we measured the resistivity but not the con- 
ductivity. However, because bismuth is a perfectly com- 
pensated semimetal, the Hall components of the conduc- 
tivity tensor are small compared to the diagonal ones. 
Therefore, the relation between the conductivity and re- 
sistivity in a strong magnetic field is very simple, namely 
pn ~ 1 / cru . Hence the oscillatory corrections are related 
by pn/pn = — crn/cii. This means that the maxima in 
(Tn, corresponding to the Fermi energy crossing the cen- 
ter of the Landau level, appear as minima in pn, as is 
observed experimentally (Fig. 1). 

In Eqn. (0) we emphasize the fact that the oscillatory 
behavior is given by the superposition of two different 
oscillations, one for holes with spin parallel to the field 
("spin up" holes) and the other for holes with spin an- 
tiparallel to the field ("spin down" holes). As has already 
been pointed out, Eqn. (Q) shows sharp features when 
the argument of one of the cosines is not p-dependent; in 



fact, in such a case v oc J^Li Vv^ ~~ * 00 • ^his ua PP ens 
when the factor multiplying 2np is an integer number. 
This singularity corresponds to the ID singularity in the 
density of states for the energy approaching the bottom 
of the Landau band. Using this observation, one read- 
ily finds the spacing between resistance minima for "spin 
down" and "spin up" holes taking the difference between 
the position of two singularities in Eqn. (Q) with opposite 
spin: this is given simply by 
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where ez is the Zeeman energy and I is an arbitrary 
integer number. The spacing in the inverse magnetic 
field is then given by 
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go is the part of the gyromagnetic ratio which contributes 
only to the fractional part of the ratio between the Zee- 
man energy and the cyclotron energy. The value of go 
can then be obtained from the fit to the quadratic de- 
pendence in Fig. (^): go — 2.68(4). The error here is 
determined by the scatter of the data. 

At this point, we need to obtain the exact value of £, 
i.e. the integer part of ez/h-^a which will be related to 
a second contribution to the gyromagnetic ratio g\ : usu- 
ally, this is done from theoretical considerations on the 
band structure but in our case it is straightforward to 
obtain this value from the observation that the last mini- 
mum in Fig. (Q) is not split. It is helpful to interpret Eqn. 
([j]) as the sum of two oscillations with the same period, 
only slightly displaced in phase due to the contribution of 
the spin, in such a way that the first peak of the doublet 
(in an increasing field) corresponds to spin down holes 
and the second to spin up holes with different Landau 
indices. The last split minimum at n = 2 corresponds to 
the crossing of the Fermi surface of the levels (v = 0, — ) 
and [y = 2,+). The spin up holes are still oscillating: 
this is the meaning of the unsplit minimum that signs the 
crossing through the Fermi surface of the level iy = 1, +): 
therefore we can conclude that the integer part of the ra- 
tio of the Zeeman energy to the cyclotron energy is equal 
to 2, which gives us g\ = 2m ai = 32.6(3). The preci- 
sion in this evaluation is primariljs-jdue to uncertainties in 
the effective mass determination. El Summing up the two 
contributions, we obtain 



.9 = 5i+.9o = 35.3(4) 



(5) 



which coincides with a ratio between the Zeeman energy 
to the cyclotron energy equal to ez/hcu c — 2J6(2). 

In agreement with previous investigationsE3, no sign of 
electronic contribution to the Shubnikov-de Haas effect 
has been identified. The electrons are expected to have a 
period three times shorter than the holes, the difference 
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being due to the non parabolicity of their energy spec- 
trum. At the same time, the temperature and the Dingle 
temperature smearing factor should have a bigger numer- 
ical factor, in such a way that the electronic oscillation 
should be depressed by a factor 3. But the only signature 
of the presence of the electrons is related to some mod- 
ulation of the hole oscillation that can be quantitatively 
interpreted assuming a Dingle temperature for electrons 
almost 10 times bigger than for the holes. On the other 
hand, no theoretical arguments can be put forward to 
explain the difference in the Dingle temperatures viz. in 
the carrier mean free time. 

In conclusion, we have shown in this work that the 
study of magneto-oscillatory phenomena at high B/T 
can reveal detailed information about electronic structure 
that is not evident at higher temperatures or lower fields. 
By high B/T we mean fields high enough and tempera- 
tures low enough such that hui c 3> fcsT. In addition the 
system must be clean enough, as are our bismuth crys- 
tals, to assure that the additional condition huj c 3> fcsTb 
or, equivalently uj c ^> 1/2-kt, holds. When both of these 
conditions are satisfied, higher harmonics in Eqn. (|l|) 
become manifest and quantitative information on the g- 
factor and spin polarization can be obtained. With in- 
creasing field the spin polarized Landau levels that cross 
the Fermi energy give rise to increasingly sharp structures 
that are limited only by the disorder broadened width 
of each Landau level. At n = 1 near a field of 9 T the 
sharpness is most pronounced (lower left insert of Fig. 1) 
with a factor of two change in resistance for a less than 
3% change in field. We have shown that at high B/T 
the appearance of these sharp features is a direct conse- 
quence of the high harmonic content of the Shubnikov-de 
Haas oscillations which together with the occurrence of 
well-defined spin-split doublets reveals important details 
in electronic structure. 
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